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ABSTRACT 

We report on XMM-Newton spectroscopy of the low-luminosity active galaxies (LLAGN) M81 and 
NGC 4579 both of which have known black hole masses and well-sampled spectral energy distributions 
(SED). The iron Ka line profiles from both the LLAGN can be described in terms of two components 
- a narrow line at 6.4 keV and a moderately broad line (FWHM ~ 2 x 10'* km s~ ) arising from highly 
ionized, He- like or H-like species {E ~ 6.8 keV). We interpret the broad lines arising from an accretion 
disk the inner edge of which is restricted to large radii (ri„ ~ lOOrg). However, the Eddington ratio, 
L/LEdd, of these sources, is 3-4 orders of magnitude lower than that required to photo- ionize a cold disk 
to He-like iron. We suggest that the lines can be explained as collisionally ionized X-ray lines arising 
from the transition region between a hot (radiatively inefficient) flow in the inner regions and a cold disk 
outside r ^ lOOrg. The accretion flow geometry probed by our XMM-Newton observations is consistent 
with the truncated disk models proposed to explain the SED of LLAGNs. 

Subject headings: accretion, accretion disks - galaxies: active ~ X-rays: galaxies 
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1. INTRODUCTION 

Active galactic nuclei (AGN) typically display a 'big 
blue bump' in their optical/UV spectra interpreted as 
blackbody emission arising from optically thick, geomet- 
rically thin accretion disks. Their hard X-ray power-law 
spectrum is thought to be produced by Comptonization 
of the soft disk photons in a hot corona above the accre- 
tion disk (e.g.; Haardt & Maraschi 1991). Low luminosity 
AGN (LLAGN), emitting well below their Eddington hmit 
{L/LEdd < 0.01), typically lack the 'big blue bump' (e.g.. 
Ho 1999; Di Matteo et al. 2003 for the case of M87) and for 
this reason it has been proposed that in these objects the 
accretion disk consists of of two zones: an outer thin disk 
that extends from some large radius down to a transition 
radius and an inner low-radiative efficiency accretion flow 
close to the black hole. 

Fluorescence Fe Ka line emission can be a very pow- 
erful probe for the accretion flow geometry around black 
holes. Such lines have been interpreted as being produced 
through X-ray irradiation of the accretion disk; there- 
fore the physical width of the line can serve as a trace 
of the inner accretion disk extent. Observations of rel- 
ativistically broadened iron Ka lines from some Seyfert 
galaxies (e.g., MCG-6-30-15 - Wilms et al 2001; Fabian 
et al. 2002), which require that the disk must be irra- 
diated from an X-ray source within 6rg, suggests that at 
least in this case the accretion disk extends to the last 
stable orbit. It is interesting to note that MCG-6-30-15 
(FWHM(Ha)~ 2200 km s"^ Sulentic et al. 1998) is hke 
narrow-line Seyfert 1 galaxies (NLSl), which are thought 
to be accreting at relatively high accretion rates. 

In contrast, Lasota et al. (1996) suggested that the ac- 
cretion disk in most LLAGN must be truncated. Gammie, 
Narayan, and Blandford (1999) have shown the broad- 
band spectrum of NGC 4258 can be explained by emission 
from a geometrically thin accretion disk + ADAF model 



with a transition radius at ~ 10 — 100 Schwarzschild radii 
(rs)- On the basis of their optical/UV continuum spectra, 
Quataert et al. (1999), have provided further evidence for 
such accretion flow geometry and showed that the spec- 
tral energy distributions of M 81 and NGC 4579 can be 
explained by an ADAF + disk with a transition radius at 
~ lOOrs. 

One prediction from any such models, as discussed in 
Quataert et al (1999), is the lack of a relativistically broad- 
ened iron Ka line in the X-ray spectrum of these objects. 
The presence, or absence, of this feature in X-ray spectra 
from Chandra and XMM-Newton can therefore be used 
to test the accretion flow geometry proposed for LLAGN. 
Specifically, any detection of relativistically broad iron Ka 
line would strongly argue against truncated disk models. 

In this paper, we utilize XMM-Newton observations of 
two of the best studied LLAGN, NGC 4579 and M 81, to 
investigate whether the observed iron Ka profiles are con- 
sistent with truncated disk models. Both NGC 4579 and 
M 81 show broad Ha line and are classified as type 1 AGN 
(Barth et al. 2001; Fihppenko & Sargent 1985). 

2. OBSERVATION AND DATA REDUCTION 

XMM-Newton observed NGC 4579 and M 81 on 12 June 
2003 and 22-23 April 2001 respectively, using he European 
Photon Imaging Camera (EPIC), the reflection grating 
spectrometer (RGS), and the optical monitor (OM) in- 
struments. The EPIC PN (Striider et al. 2001) and MOS 
(Turner et al. 2001) cameras were operated in full frame 
mode using the thin filter for the observation of NGC 4579. 
For the observation of M 81, the EPIC PN was operated 
in small window, the MOSl camera was operated in fast 
uncompressed mode and the M0S2 camera was operated 
in full frame mode, with the medium filter in place for 
all three detectors. We present the spectral data obtained 
from the EPIC cameras for NGC 4579 and the spectral 
data from the EPIC PN and M0S2 camera for M 81. The 
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EPIC PN exposures were 22ks and 130ks for NGC 4579 
and M 81, respectively. 

The raw events were processed and filtered using the 
most recent updated calibration database and analysis 
software (SAS v5. 4. 1) available in August 2003. Events in 
the bad pixels file and those adjacent pixels were discarded. 
Only events with pattern — 4 (single and double) for the 
PN and — 12 for the MOS were selected. Light curves ex- 
tracted from source free regions were examined for flaring 
particle background and affected exposure intervals were 
excluded from the rest of the analysis. This resulted in 
'good' exposure time of 15.2 ks and ks for NGC 4579 and 
71.9 ks M 81, respectively. 

Source spectra were extracted from a circular region of 
radius 50" for both NGC 4579 and M 81. Associated back- 
ground spectra were extracted from source free regions. 
The total number of source photons, detected with the PN 
is 4.1 X 10^ for NGC 4579 and 5.44 x 10'^ for M 81. The 
source spectra were grouped to a minimum of 20 counts 
per spectral channel with appropriate sampling of the data 
(~ 3 channels per spectral resolution elements) and were 
analyzed using XSPEC 11.2.0 (Arnaud 1996). All the er- 
rors quoted below were calculated at 90% confidence level 
for one interesting parameter. 

3. DETECTION OF COMPLEX IRON LINE PROFILES 

In the following analysis, only the EPIC data above 
2.5 keV are examined, the results presented in this paper 

are not critical to the soft X-ray spectrum. The soft X-ray 
spectrum of M 81 obtained from the RGS observation has 
been studied by Page et al. (2003b). For each object, at 
first the MOS and PN spectra were fitted separately to 
check for the uncertainties due to cross-calibration prob- 
lems. Wc found generally good agreement between M0S2 
and PN cameras (AF 0.03) in the 2.5 - 10 keV bands. 
However, the MOSl data results in a flatter (AF ~ 0.15) 
spectrum. These discrepancies are consistent with that 
found by Molendi & Sembay (2003). Therefore we present 
the spectral results obtained by fitting the same model 
jointly to the PN and M0S2 data in the 2.5 - 12 keV 
band, and allowing the relative normalizations to vary. 

3.1. NGC 4579 

We fit the PN and M0S2 spectrum of NGC 4579 jointly 

with a power-law model modified by photoelectric ab- 
sorption. This fit resulted in a of 296.6 for 253 de- 
grees of freedom (dof). The best-fit absorption column 
density is consistent with the Galactic column density 
{Nh = 3. 1 X 10^° cm-2 along the line of sight of NGC 4579. 
From this simple fit, at least two iron line components are 
evident in the X-ray spectrum of NGC 4579. Figure 1 
shows the iron Ka line profile of NGC 4579 as the ra- 
tio of the observed PN data to the the best-fit continuum 
model (see below). Addition of a narrow Gaussian line at 
^ 6.4 keV improved the fit significantly (x^ = 267.1 for 
251 dof). However, excess emission is seen at ~ 6.7 keV. 
The addition of a broad Gaussian line (model 1: power-law 
-|- narrow Gaussian + broad Gaussian) improves the joint 
fit significantly (x^ = 242.5 for 249 dof), over the best-fit 
single Gaussian line model and also results in a better fit 
than the best-fit from a model consisting of an absorbed 
power-law and two narrow Gaussian lines (x^ = 257.0 for 



249 dof). The best-fit parameters are F = 1.77to o7, nar- 
row Gaussian: EKa,n = 6.40 keV (fixed), EW = 106 eV, 
broad Gaussian: EKa,b = 6.79t;^i3 keV, a = 231^9^ eV, 
EW = 287 eV (see model 1 in Table 1). 

The observed flux of NGC 4579 in the 2.5-12 keV band 
is 3.9 X 10~^^ erg cm~^ s~^, and the corresponding lumi- 
nosity is 1.2 X 10*° erg s~^ (for a distance of d = 16.7 Mpc 
as measured by Pierce & TuUy, 1998). An ASCA ob- 
servation of NGC 4579 in 1995 showed a strong (EW 
~ 500 eV), marginally broad, symmetrical, iron Ka line 
profile at 6.73 ± 0.13 keV, while the previous 1998 ob- 
servation showed a narrow line at 6.39 ± 0.09 keV (EW 
^ 260 eV) with no clear line at ^ 6.7 kcV. It appears 
that both the above components observed with ASCA are 
clearly detected in the current XMM-Newton observation. 

3.2. M 81 

The PN and M0S2 spectra of M81 were also fit by 
a power-law modified by photoelectric absorption. The 
best-fit absorption column density is consistent with the 
Galactic value (TVh = 4.2 x 10^° cm'^). This simple 
fit is statistically imacceptable (x2 of 397.5 for 303 dof) 
due to significant spectral features in the iron Ka region 
(see Fig. 1). The complex line profile of M 81 is not 
well described by a narrow Gaussian line with the power- 
law plus narrow Gaussian model resulting in a poor fit 
(x^ = 363.0 for 303 dof). The addition of a second narrow 
Gaussian line at ~ 6.7 keV improves the fit significantly 
(X^ = 318.3 for 301 dof). Varying the width of this line 
further improves the fit (x^ = 299.3 for 300 dof). The 
best-fit parameters are Fx = 1.87±0.02, narrow Gaussian: 
EKa,n = 6.40 keV (fixed), EW = 31 eV, broad Gaussian: 
EKc^,b = 6.79tg-2? keV, a = 188,tlr eV, EW = 101 eV. 

The observed flux of M 81 in the 2.5 - 12 keV band is 
1.2 X 10""'^^ erg cm^^ and the hmiinosity in the same 
band is 1.8 x 10^" erg s^^ (for a distance of d = 3.6 Mpc; 
Frcedman ct al. 1994). 

The XMM-Newton observations of M 81 confirm the 
previous ASCA and BeppoSAX results that detected an 
intrinsically broad or complex iron line centered at ~ 
6.7 keV, with an equivalent width of ~ 170 ± 60 eV 
(Ishisaki et al.l996). Neither ASCA nor XMM-Newton 
observations detected an absorption edge at ^ 8.6 keV, 
similar to that observed previously with BeppoSAX (Pel- 
legrini et al. 2000). 

4. INFERENCES FROM THE IRON LINES 

Both NGC 4579 and M 81 show iron line emission with 

broad {FWHM ^ 20000 km s""'^), symmetrical line pro- 
files, in addition to narrow iron Ka lines at 6.4 keV. 

The broad iron Ka lines cannot be easily interpreted 
as thermal emission from the host galaxy because the in- 
ferred thermal broadening from such a system is too small 
(FWHM ^ 350 km s~^ for a temperature of 10^ K) to 
explain the widths of the observed lines. A thermal origin 
for the lines is further ruled out by the variability of the 
lines as observed by ASCA (Terashima et al. 2000). 

The observed broad line from NGC 4579 and M 81 may 
thus arise from the accretion disk, as generally assumed in 
Seyfert galaxies. The broad Gaussian lines in the current 
models were replaced with the disk-line model of Fabian et 
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al. (1989) in order to test this hypothesis. This model as- 
sumes a Schwarzschild geometry and a disk cmissivity that 
is a power-law as a function of disk radius (e oc r^'^). The 
emissivity index was fixed at g = 2.5, the average value for 
Scyfcrt 1 galaxies derived by Nandra ct al. (1997). The 
outer disk radius was also fixed at lOOOrg and the disk in- 
clination angle i was fixed at 30 deg. The free parameters 
of the disk-line model were the line energy, the inner radius 
of the accretion disk, and the line flux. This model (model 
2: power-law + narrow Gaussian + disklinc) results in 
an acceptable fit for both AGN; x^/dof = 240.6/249 for 
NGC 4579 and x^/dof = 296.3/300 for M 81. The contin- 
uum and narrow Gaussian line parameters did not change 
significantly (see Table 1). 

Figure 2 shows the observed PN data, the best-fit model 
(which includes an absorbed power law a narrow Gaus- 
sian and a diskline) and the deviations (in units of a) for 
NGC 4579 and M 81. The iron Ka profiles, shown in Fig. 1 
as the ratio of observed data and the best-fit continuum 
model, were obtained by setting the normalizations of the 
narrow and broad Ka lines to be zero. Only the ratios to 
the PN data are plotted in Fig. 1 for clarity. In order to 
evaluate the dependence of the line properties on disk in- 
clination, we also fit the data for i = 60 deg. This results in 
Ediskiine = 6.70^^9 keV; nn = 232t*?°rg for NGC 4579 
and Ed^skHne = 6.821^0 kcV; r,„ > leOr, for M 81. Fit- 
ting accretion disk models to the data suggest that the line 
emitting region of the disk is truncated at an inner radius 
of about or greater than lOOfg (the narrow core to the 6.4 
keV line suggests that this component must also originate 
at a substantial distance from the black hole). 

Our best-fit model does not show any residuals indi- 
cating the presence of a substantial reflection continuum 
(from cither neutral of partially ionized matter from a disk 
close to the black hole) in the spectral data of either M81 
or NGC 4579 (sec Fig. 2). The apparently lack of reflec- 
tion continuum was also confirmed by the detailed spectral 
analysis of the BeppoSAX observation of M81 (Pellegrini 
2000). It is therefore imlikely the the 6.7 keV emission 
is produced by reflection from an ionized disk (see also 
§5.2). To assess the possible origin of the ^6.7 keV broad 
feature in the spectrum of NGC 4579 as the Compton re- 
flection emission, we fltted a Compton reflection model 
(pexrav; Magdziarz & Zdziarski 1995) along with the ab- 
sorbed power-law and narrow Gaussian line model to the 
PN spectrum of NGC 4579. The best-fit of this refiection 
model results in an unrealistically large relative reflection 
parameter, R = 8.2l^ 5 (where R = 1 corresponds to a 
solid angle of 2tt subtended by the reflector) and the fit 
worsens (Ax^ = -1-20.7 for an additional dof) in compari- 
son to the best-fit obtained with the absorbed power law 
and double Gaussian line model (see Table 1). Thus the 
unphysical value of the refiection parameter and the poor 
quality of the fit rule out Compton refiection emission in 
place of the fluorescence iron Ka emission. 

Finally, a combination of multiple narrow Gaussian lines 
was also tested as a model to fit the complex iron line 
profiles observed in these sources. Three narrow Gaus- 
sian emission lines were added to the absorbed power-law 
model with fixed line energies of 6.4, 6.7 and 7.0 keV, cor- 
responding to cold. He-line, and H-like iron. The fits based 
on this model are statistically acceptable fits (NGC 4579, 



= 248.5 for 250 dof; M 81, x^ = 304.1 for 301 dof). 
These are worse fits compared to the diskline models 
(Ax^ = +7.9 and -1-7.8 for NGC 4579 and M 81, respec- 
tively) at a confidence level of 99.5% based on an F-test. 
With the current generation of X-ray spectroscopy satel- 
lites it is not possible to clearly distinguish between the 
multiple (> 3) narrow line and the disk-line models how- 
ever, observations with. Astro E2, and Constellation- X, 
with their correspondingly higher spectral resolution and 
sensitivity will be particularly important to further inves- 
tigate the complexity of these iron line profiles. 

5. ORIGIN OF IRON LINE COMPLEX 

Both NGC 4579 and M 81 show complex Iron Ka fea- 
tures which are well described by a combination of a nar- 
row Gaussian and either a broad Gaussian emission line, a 
disk-line, or a combination of narrow emission lines from 
multiple ionization states of iron. 

5.1. The neutral, narrow iron Ka line 

Despite the good spectral resolution of the EPIC cam- 
eras, the line energies and the widths of the observed nar- 
row Ka lines are not well constrained. The line parameters 
for NGC 4579 are similar to the narrow lines commonly 
observed from Seyfert galaxies with Chandra and XMM- 
Newton (e.g., Yaqoob et al. 2001; Pounds et al. 2001), 
however, the strength of the narrow Ka line is weaker in 
M 81, in comparison to that of 'typical' Seyfert 1 galaxies 
(see e.g. Page et al. 2003a). The observed narrow-lines 
from NGC 4579 and M 81 can be produced anywhere in the 
broad-line region, dusty torus, outer regions of an standard 
accretion disk or a disk with a flat emissivity law (Yaqoob 
et al. 2001,2003). The weaker narrow fine of M 81 may 
suggest a smaller covering fraction of the material respon- 
sible for the narrow line emission in NGC 4579 than that 
in M 81. 

5.2. The ionized, broad iron Ka line 

The widths and line center energies of the broad Ka 

lines from NGC 4579 and M 81 are remarkably similar 
{E - 6.7 - 6.8 keV; FWHM ~ 20000 km s"^) suggestive 
of both an origin in similar regions the AGN and a sim- 
ilarity in the ionization stages (He- or H-like iron) of the 
emitting material. The observed line energies are higher 
than those generally observed from Seyfert 1. In general, 
only narrow-line Seyfert 1 galaxies (NLSls), with steep 
X-ray spectra, regularly show strong, highly ionized iron 
Ka lines (Dewangan 2002). In such cases the lines are 
interpreted as being produced by reflection from mate- 
rial in the disk photo-ionized by the strong X-ray con- 
tinuum. The physical problem with explaining the lines 
in M81 and NGC 4579 in a similar way is that we do 
not expect photoionization to be important in these low 
luminosity sources. For photo-ionized material, the ion- 
ization parameter can be written as ^ = where L is 
the luminosity of ionizing photons, n is the density of the 
material, and R is the distance of the material from the 
ionizing source (Kallman & McCray 1982). For standard 
thin disks. Matt et al. (1993) have shown that the ioniza- 
tion parameter is strongly dependent on the accretion rate 
relative to the Eddington rate and hence can be rewritten 
as, ^ oc (L/LEdd)^- In order to ionize He- like iron, ^ should 
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be at least ~ 500 or L/LEdd > 0.2 (Matt et al. 1993). This 
result is consistent with the observation of highly ionized 
iron Ka lines from some NLSls as they are thought to emit 
a higher fraction of their Eddington luminosity (Pounds, 
Done, & Osborne 1995). For NGC 4579, the total nuclear 
luminosity of ~ 5 x 10"^^ erg s~^ (Barth et al. 1996) and 
black hole mass in the range of 1.8 - 9.3 x lO'^ M0 (Barth 
et al. 2001) correspond to L/LEdd ~ (0.5 - 2.2) x 10"* 
while Laor (2003) reported L/Lsdd ^ 1-35 x lO"''. Ho et 
al. (1996) found L/LEdd (2 - 10) x 10"^ for M 81 which 
has a black hole mass of 4 x IQ^Mq. Thus the ionizing 
luminosity of both the LLAGNs NGC 4579 and M 81 are 
3-4 orders of magnitude smaller than that required for a 
photo-ionized accretion disk with He-like iron. Clearly, the 
dominant ionization process in the disk of NGC 4579 and 
M 81 cannot be photo-ionization. 

The widths of the ionized iron Ka lines from NGC 4579 
and M 81 (FWHM- 20000 km s"^) arc not as broad as 
the line widths predicted by models of an accretion disk ex- 
tending down to the last stable orbit. Assuming Keplerian 
motion of the emitting material (i.e., r ~ (c^/^^FVFi?M)^s) 
the size of the broad iron Ka line emitting material is es- 
timated to be - 150rg for NGC 4579 and - 230rg for 
M 81, independent of black hole mass. Thus the majority 
of iron Ka photons are not produced near the last stable 
orbit for a Schwarzschild (rj„ = 6rg) or a Kerr black hole 
(1.235rg < Tin 6rg). The same conclusion can be drawn 
from the disk-lino fits which resulted in rj„ = 50i2i^''g for 
NGC 4579 and 65t2?^rg for M 81. 

Two possible scenarios that can explain the observed 
moderate line widths and their ~ 6.7 — 6.8 keV centroid 
energy are: {i) the accretion disk is fully ionized for r < r^^ 
so there are no lines arising from within this region, or {ii) 
the accretion disk is truncated at r ^ ri„. As discussed 
above, however, the accretion disk cannot be photo-ionized 
by irradiation, and thus must be intrinsically hot. Both 
(i) and {ii) then imply that the inner regions of the accre- 
tion flows in M 81 and NGC 4579 LLAGNs are likely to 
be different from those of typical Seyfert galaxies. 

5.3. The truncated disk model 

Our results are consistent with the suggestion that 
LLAGN have accretion disks that are truncated at some 
large radii and have a hot flow in the inner regions (La- 
sota ct al. 1996; Narayan 1996; Gammic et al.; Li & 
Wang 2000). In particular, the inner radius Tin ^ lOOr^ of 
the accretion disk implied by the observed iron line width 
is fully consistent with the position of the truncation ra- 
dius deduced from the modeling of the optical/UV spec- 
tra of M 81 and NGC 4579 by Quatacrt ct al. (1999; 
assuming a geometrically-thin accretion disk and an opti- 
cally thin, two-temperature ADAF in the inner regions). 
Note also that in this model, the continuum X-ray emis- 
sion is produced in the inner region of the hot flow (from 
r — 1 — lOOrj;) and is due to Comptonization of synchrotron 
emission from the flow. 

We suggest that the relatively broad 6.7 keV lines may 
arise from regions of the flow close to or at the transi- 
tion radius, where the hot, fully ionized flow starts cool- 
ing to connect to the outer optically thick, geometrically 
thin region. Using the dynamical solutions of the tran- 
sition regions from a thin disk to an ADAF recently de- 



rived by Manmoto et al. (2001) we can try to estimate 
the densities and temperatures typical for this transition 
region. In order for the hot flow to connect to a stan- 
dard thin disk the temperature and velocity of the gas 
must change by orders of magnitude over a small transi- 
tion region between the ADAF and the thin disk. For a 
transition radius of Vtr ^ lOOr^, the temperature (T oc c^) 
scales as ~ in the inner regions, and as it approaches 
the transition radius, it then drops sharply from values of 
T - lO'^ K at r - 50r, to - 10*^ K at r - lOOr,, and 
cooler beyond. As pressure balance has to be maintained, 
the density of the flow increases sharply (by ~2 orders 
of magnitude) in the same range of radii, reaching values 
of rie ^ 2 X 10^^ cm~^ at r lOOr,,. Even using the 
values of the flow density derived above, photo-ionization 
in the transition region is not likely to be important be- 
cause of the low level of the ionizing hard X-ray flux, in 
fact, at the transition region wo find that the ionization 
parameter ^ ~1. Instead, the observed He-like and H- 
like iron lines are likely to represent coUisionally ionized 
X-ray lines produced in a hot, T ~ 10^ K, transition re- 
gion. Specific calculations of the emission line profiles from 
this region are not available at the moment and would be 
necessary to fully explore this hypothesis. To test this 
inference, therefore, we consider a simple thermal plasma 
model (mekal; Liedahl, Osterheld, & Goldstein (1995) and 
references therein) and fit it in combination with a power- 
law continuum and a narrow Gaussian omission line at 
6.4 keV (fixed) to the PN and MOD spectra jointly. The 
density of the plasma was fixed to n = 2 x 10^^ cm~^ 
(as derived above) and no intrinsic absorption is required. 
The mekal model was convolved with a Gaussian in or- 
der to mimic the lino broadening due to Keplerian mo- 
tion of the plasma at the transition region. The Gaussian 
fT was fixed at 0.2 keV, similar to that derived for the 
broad Gaussian component in the double Gaussian model 
for the iron line (see Table 1). The model provided a 
good fit to the data {x^/dof = 246.4/249 for NGC 4579 
and 298.0/296 for M 81), albeit somewhat worse than the 
models incorporating either a broad Gaussian or diskline. 
As an example, this model is shown in Figure 3 for the PN 
data of M 81. The best-fit parameters are F = 1.73t2:??, 
kT = S.StIi keV for NGC 4579 and T = 1.87 ± 0.04, 
kT = 8.5^3 keV for M 81. The best-fit plasma tem- 
peratures for both the AGN are similar to that deduced 
above for the transition region. The power-law and the 
mekal fluxes are 2.6 x 10"-^^ and 1.3 x 10"-^^ erg cm~^ 
for NGC 4579 and 1.0 x 10"" erg cm-^ and 1.3 x 
10~^^ erg cm~^ s~^, respectively. 

The success of this model in reproducing the observed 
iron Ka profiles suggests that the observed lines could be 
produced by collisional excitation in a relatively hot tran- 
sition region between a hot flow and a thin disk. The 
truncation radii of ~ lOOrs radii inferred by Quatacrt et 
al. (1999) to fit the optical/UV to X-ray SED is some- 
what larger than the best-fit inner radii inferred from 
the diskline model fits to the iron Ka line profile. How- 
ever, is has been noted the if a standard disk is con- 
nected to an ADAF across a narrow transition layer, the 
gas near the transition radius must rotate with a super- 
Keplerian angular velocity (Abramowicz, Igumenshchev, 
& Lasota 1998) This means vfwhm(Pg Ka) > c{rg/rtr) 



or rtr > Tin ~ {'^'^ /^FWHM)^g- Additionally, the transi- 
tion region, being hot (~ 10^ K), may not contribute the 
optical/UV emission which determined the truncation ra- 
dius quoted by Quataert et al. (1999). Note also that 
ADAFs arc likely to develop strong outflows (Blandford 
& Begelman 1999) or convection (Quataert & Gruzinov 
1999; Narayan et al. 2000) which are likely to change the 
dynamics of the transition region. 

5.4. Alternative scenario - multiple narrow lines? 

The broad iron line components of NGC 4579 and M 81 
can also be described in terms of multiple narrow emis- 
sion lines from ionized iron, although with a lesser statis- 
tical significance. Such lines can arise from (i) hot ther- 
mal plasma, if present, outside the accretion disk e.g host 
galaxy, or (ii) a photo-ionized accretion disk with flat emis- 
sivity law extending down to the last stable orbit (Yaqoob 
et al. 2003). The former process has already been ruled 
out due to the variability of the iron line from NGC 4579 
(Terashima et al. 2000). The later possibility is unlikely 
due to a low ionization parameter (see §5.3). Thus the 
case for multiple narrow lines for the iron Ka complexes 
of NGC 4579 and M 81 is poor at present. 

6. CONCLUSIONS 

Using the XMM-Newton EPIC observations, we de- 
tected complex iron line profiles from two LLAGN 
NGC 4579 and M 81. The line profiles are well described 
by a combination of a narrow neutral line and a broad 
Gaussian (FWHM ^ 20000 km s~^) or an accretion-disk 
line arising from highly ionized (He- or H-like) iron. The 
line profiles from both the AGN are remarkably similar 
except that the neutral iron Ka line from M 81 is much 
weaker than that that from NGC 4579. The main result 
of this work is that the broad and highly ionized iron Ka 
lines from two LLAGNs NGC 4579 and M 81 are consis- 
tent with accretion disks whose inner edges are restricted 
to large radii ri„ ~ lOOr^. 

We thank an anonymous referee for the comments and 
suggestions on this paper. This work is based on obser- 
vations obtained with XMM-Newton, an ESA science mis- 
sion with instruments and contributions directly funded by 
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Fig. 1. — Iron Ka line profiles of NGC 4579 and M 81 shown as the ratio of the observed PN data and the best-fit continuum model (see 
text). 
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Fig. 2. — The observed EPIC PN and MOS2 spectra of NGC 4579 and M 81, the best-fit models consisting of an absorbed power-law, a 
narrow Gaussian, and a diskline (see text). Lower panels show the deviations of the observed data from the best-fit models. 




Fig. 3. — The best-fit model consisting of (i) mekal plasma model, (ii) a narrow Gaussian, and (iii) a power-law derived from the joint fit 
to the PN and MOS2 spectra of M 81. The model is shown for the PN data only for clarity, the model for the MOS2 data differs in relative 
normalization by a small fraction (~ 0.02%). The mekal model was convolved with a Gaussian of cr ~ 0.2 keV (see text) 
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Table 1 

Best-fit spectral model iwii ameters for NGC 4579 and M 81. 



Comp. 


Parameter'' 


NGC 4579 


M 81 






model 1*" 


model 2^ 


model l'' 


model 2^ 


PL 


Tx 


1 77+U.U6 
' ' -0.07 


1 77+0.U6 
'-■ ' ' -0.07 


1.861'^:^^ 








J--UZ_Q Qg 


1 02+°" 




3.612:1 


Gaussian 


EkcA keV) 


6.40(f) 


6.40(f) 


6.40(f) 


6.40(f) 




a ( keV) 


0.00(f) 


0.00(f) 


0.00(f) 


0.00(f) 




f 


A 7+2-3 


4 

^••J-1.9 


7+1.6 

■J-'-i.g 






EW ( cV) 


106 


96 


31 


35 


Gaussian 


EkcA keV) 


6.79l°:li 




6.79«:°? 






a ( eV) 


23llif 




188il« 






f 


Q Q+4-4 
^•^-3.8 










EW ( eV) 


287 




101 




DL 


Edisklinei keV) 




f! 70+0. 06 
0- '8-0.09 




6.84ir;.?,^ 




q 




2.5 




2.5(f) 




n/rg 








104^^1" 




Vo/rg 




1000(f) 




1000(f) 




i 




30(f) 




30(f) 








10.01?,:^ 




10.ll|°3 




EWd{ eV) 

Xmin /dof 




302 




104 


Total 


242.5/249 


240.6/249 


299.3/300 


296.3/300 



"(f) indicates a fixed parameter value. 

''model 1: absorbed power law + narrow Gaussian + broad Gaussian; model 2; absorbed power law + narrow Gaussian + diskline. 
"=10-3 photons cm-2 s"! keV"! at 1 keV. 
■^lO-"^ photons cm-2 s'^ 



